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Particle  Size  Distributions  during  Diffusion  Controlled  Growth  and  Coarsening 

Oleg  N.  Senkov 

UES,  Inc.  4401  Dayton-Xenia  Rd.,  Dayton,  Ohio  45432 

ABSTRACT 

Theoretical  analysis  of  particle  growth  and  coarsening  conducted  in  the  present  work,  which 
includes  stochastic  flux  of  a  solute  inside/from  the  growing  particles,  predicts  normal  distribution 
of  the  particle  sizes,  independent  on  the  particle  volume  fraction.  This  result  agrees  well  with  the 
numerous  experimental  data  and  it  solves  a  mysterious  problem  related  with  a  non-symmetrical 
size  distribution  in  earlier  Lifshitz-Slyozov-Wagner  theory. 

INTRODUCTION 

Growth  and  coarsening  of  second  phase  particles  in  a  supersaturated  solution  are  important  parts 
of  phase  transformations  which  tune  various  physical,  chemical  and  mechanical  properties  of 
materials.  These  transformation  processes  occur  when  the  nucleation  stage  of  the  second  phase 
particles  has  completed  and  further  decrease  in  the  free  energy  of  the  system  is  mainly  due  to  a 
decrease  in  the  total  interfacial  energy.  The  particle  growth  originates  from  supersaturation, 
which  is  still  present  in  the  system  after  completion  of  the  nucleation  stage,  and  from  the 
concentration  gradients  around  the  particles  of  different  sizes  caused  by  the  capillary  effect  (the 
Gibbs-Thomson  equation),  which  result  in  preferential  growth  of  larger  particles  and  shrinkage 
of  smaller  particles  and  in  a  continuous  increase  in  the  particle  average  size  [1,2].  The  major 
advance  for  the  theoretical  description  of  this  process  was  made  by  Lifshitz,  Slyozov  and 
Wagner  (so  called  LSW  theory)  [3,4].  This  theory  describes  very  well  evolution  of  averaged 
growth  parameters,  such  as  the  average  radius,  r  ,  and  the  number  density  of  the  particles,  N  , 
with  time  t.  For  example,  for  the  diffusion  controlled  growth  process,  the  LSW  theory  predicts 


the  following  behavior: 

r  3  -  L03  =  Xlsw  t 

(1) 

N  1  -  N0 1  =  (4tc/3/v)  Xlsw  t 

(2) 

1 


where  r0  and  N0  are  the  average  radius  and  the  number  density  of  particles  at  time  t  =  0,  /y  is 
the  particle  volume  fraction  and  Alsw  is  the  rate  constant,  which  is  given  by: 


^LSW  - 


W„,c 

9  RT 


(3) 


In  Equation  (3),  y  is  the  particle-matrix  interface  energy,  D  is  the  diffusion  coefficient,  Vm  is  the 


molar  volume  of  the  second  phase  and  cK  is  the  equilibrium  saturated  concentration  of  the 


alloying  element  in  the  matrix  at  the  temperature  T,  and  R  is  the  gas  constant.  Many  experimental 
results  on  the  particle  growth  and  coarsening  indeed  follow  these  relationships.  However,  the 
LSW  theory  leads  to  a  highly  asymmetric  particle  size  distribution  (PSD): 


0LSW  : 


3  +  pj  U-2  p 


exp  -  at  r  <  1.5  and  (Plsw  =  0  at  p  >  1.5 

13-2 p) 


with  a  sharp  cutoff  near  the  particle  radius  r=  1 .5  r  ,  which  virtually  does  not  agree  with  any 
experimentally  reported  particle  size  distributions  [2,5].  The  reported  distributions  are  generally 
broader  and  more  symmetric  than  the  LSW  prediction  (Figure  1). 

A  major  problem  with  the  LSW  theory  is  believed  [1,2]  to  be  due  to  using  the  mean  field 
approximation  for  the  kinetic  equation.  This  approximation  assumes  that  the  growth  of  a  particle 
does  not  depend  on  diffusional  interactions  between  nearest  particles  and  it  does  not  take  into 
account  the  effect  of  the  volume  fraction  of  particles  on  the  rate  of  growth  of  these  particles. 

After  taking  into  account  the  volume  fraction  of  the  particles  and  their  soft  interactions,  modified 
LSW  theories  [6, 7, 8, 9,  see  also  ref.  2  for  review]  made  the  rate  constant  in  Equation  (1) 
dependent  on  the  particle  volume  fraction  and  improved  PSDs  by  making  them  broader  and  more 
symmetric.  However,  these  modified  PSDs  become  near  symmetrical  only  at  high  particle 
volume  fractions  (fy  >  0. 1)  [2],  The  fact  that  experimental  PSDs  are  also  symmetric  at  very  low 
fv  may  indicate  that  some  other  quantity  or  condition  of  the  system,  which  has  not  been  taken 
into  account  in  the  particle  growth  theories,  is  responsible  for  the  near-symmetric  distribution.  In 
this  paper  we  will  show  that  the  missed  link  between  experimental  and  theoretical  PSDs  was 
stochastic  solute  diffusion  that  was  not  taken  into  account  in  those  theories.  To  derive  the  PSD 
for  the  growth  and  coarsening  process,  we  use  a  classic  approach  developed  by  Einstein  more 
than  100  years  ago  when  he  solved  the  problem  of  Brownian  motion  [10]. 
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ANALYSIS 

Here  we  consider  diffusion-controlled  growth  and  coarsening  of  second  phase  spherical  particles 
in  a  supersaturated  solution,  although  similar  approach  can  also  be  applied  to  the  reaction-limited 
growth.  The  particles  are  distributed  by  theirs  sizes  (radii)  and  their  number  density  is  N,  which 
are  time  dependent.  The  initial  concentration  of  the  solute  element  in  the  solution  is  c0,  current 
concentration  is  c,  and  its  concentration  in  particles  is  assumed  to  be  constant  and  equal  cp.  In 
general  case,  growth  of  these  particles  is  due  to  (a)  the  diffusion  current  of  solute  across  the 
particle/matrix  interface  caused  by  the  solute  concentration  gradients  near  the  particles  and  (b) 
stochastic  diffusion  of  solute  across  the  interface  through  the  Brownian  motion.  The  first  process 
was  well  described  in  the  previous  theories,  while  the  second  process  was  ignored.  If  this 
stochastic  process  is  taken  into  account,  the  diffusion  current,  J,  of  solute  across  the 
particle/matrix  interface  per  unit  area  and  the  rate  of  particle  growth,  dr/dt,  are  given: 


J=  D 


dc 

dR 


R=r 


+  g(r,t)c(t) 


dr 

dt 


rc  pkT 


+  g(r,t)lcp% 


(5) 

(6) 


Here  R  is  the  distance  from  the  center  of  a  particle  of  radius  r,  y  is  matrix/particle  interface 
energy,  Q  is  the  particle  atomic  (or  molar)  volume,  rc  is  a  critical  radius,  given  by  a  linearized 


Gibbs-Thompson  equation,  c  =  ca 


1  + 


2]0, 

rkT 


at  which  a  particle  of  radius  rc  is  in  equilibrium 


with  the  solution  of  concentration  c{t),  c-f!  is  the  equilibrium  solute  concentration  (at  t— >oo),  t  is 
time,  k  is  the  Boltzmann  constant  (or  universal  gas  constant),  T  is  the  absolute  temperature, 
g(r,t)£;  is  a  Langevin  (stochastic)  flux,  and  E,  is  the  rapidly  fluctuating  random  term,  which  has 
the  following  properties  [11,12]: 


<%$)>  =  0  and  <^MO>  =  5 (t-f) 


(V) 
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Conditions  (7)  indicate  that  the  average  stochastic  flux  is  equal  zero,  while  its  correlation 
function  is  non-zero  and  all  the  fluctuations  at  different  times  are  independent.  In  equation  (7), 

<  >  means  the  mean  value  and  8  is  a  delta  function.  Because  of  the  stochastic  nature  of  g(r,t)^(t), 
this  value  should  not  affect  the  time  dependences  of  the  average  values  of  radius,  number  density 
and  volume  fraction  of  the  growing  particles.  However,  the  presence  of  the  stochastic  flux  will 
change  the  particle  size  distribution,  as  it  is  shown  below. 

In  the  presence  of  the  diffusion  current,  Equation  (5),  the  radius  r  of  an  individual  particle  during 
a  very  short  period  of  time,  d r,  will  increase  by  the  amount  of  A,  which  can  be  defined  from 
Equation  (6)  as: 

2iO.Dc  f  1  l'] 

A  = - —=r - dr  +  g(r,t)/cpdW(t)  (8) 

rcpkT  {  r  r) 

t 

where  W(t)  =  j  £(t)dt  is  the  Wiener  process,  which,  in  Ito  calculus  for  stochastic  processes,  has 

o 

the  following  properties:  drdW(t)  =  0  and  [dW(t)]~  =  <7  r  [  1 2],  For  each  particle,  A  has  a  different 
(positive  or  negative)  value  and  the  number  of  the  particles,  diV,  which  radius  escalations  are 
between  A  and  A+dA,  is  given  by 

dN=N</{A)dA  (9) 

where  the  function  <p{A)  is  only  different  from  zero  for  very  small  values  of  A  and  satisfies  the 

oo 

condition  J  </>(A)dA  =  1 .  If  only  stochastic  flux,  g(r,t)£,  were  present,  then  (/A  A)  =  (/A- A). 

-oo 

However,  the  presence  of  a  directional  flux  due  to  the  concentration  gradient  (Equation  5)  is  the 
main  contributor  to  the  particle  growth  and,  therefore,  (fAA)  ^  (/A- A). 

We  now  define  a  distribution  function  as  a  number  of  particles  per  unit  volume  of  radius  r 
at  time  t.  Equation  (9)  allows  us  to  determine  the  distribution  of  particles  at  time  t  +  dr  from  the 
distribution  at  time  t.  Indeed,  the  number  of  particles  at  time  t  +  dr  that  have  radii  between  r  and 
r  +  dr  is: 

oo 

d{r,t+df)dr  =  dr  J  ®(r  -  A,t)<f)(A)dA  (10) 

-oo 

Since  dr  and  A  are  very  small,  we  can  set  Taylor  series: 
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( tfr,t+z)  =  (P(r,t)  +  (dr) 


50  (r,t) 
dt 


(id 


=  Hr.l)-  A^l  +  (12) 

dr  2  8r2 

We  left  the  third  term  in  Equation  (10)  because  A  contains  the  Wiener  process  (see  Equation  8) 

2  2 

and,  therefore,  A  =  (g/cp)  (dr).  We  can  use  Equation  12  under  the  integral,  because  only  small 


values  of  A  contribute  to  Equation  10: 
<S>(r,t)+(dT)  =d\r,t)  j (/)(A)dA 


50  (r,t) 
dr 


J  A(j)(A)dA 


+ 


1  52O(r,0 

2  dr 2 


J  A2<f>(A)dA  (13) 


Taking  into  account  that  J  (f)(A)dA  =  1  and  setting 


lim  —  f A</>(A)dA  =  —  and  lim  —  f A2 </>(A)dA  = 

dz-^> 0  A  t  J  At  dz-^ti  A  t  J 


d  z 


dr 

dt 


1 


<**■-»<>  dr 


da ' 
dt 


(14) 


dr  2  .....  . 

where  —  is  the  average  particle  growth  rate  and  a"  is  the  size  distribution  variance,  we  obtain 
dt 


the  differential  equation  for  the  PSD  function  0(r,/): 


50(r,t) 


dr  50(r,t)  1  da  5“0(r,t) 


dt  dt  dr  2  dt 

This  equation  has  solution: 

f 

O  (r,t)  =  A  ex p 


dr" 


(15) 


(r  -  r  f 


A  = 


V 

2  N(t) 


2a' 


(16) 


A 


na' 


1  +  erf 


f  -  W 
r 


V  V2 a~  J 


GU 

where  A  is  a  normalizing  parameter  calculated  under  the  condition  that  J  O (r,t)dr  =  N(t). 


_  #  2 

The  values  r  and  a  in  Equation  (16)  can  now  be  correlated  with  <r>  and  <r~>\ 

|~2~  crexp(-r2/2cr2) 
n  1  +  erf  {f / CT1/2) 


<r>=  r  + 
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<r2>  =  r 2  +  cr2  +  J — 


IY  rcrexp(-r2/2cr2) 
v  7T  1  +  erf  (r/ a  yY'j 

2  crexp(-F2/2cr2)  „  „  _  ,  IY  raexp(-r2/2a2) 

- v  /  '/<  0.007 r  and  J - n/  f 


At  dr  <  0.4,  . I - v  /_  — — 

V  n  1  +  erf[r/<jf2 


tc  1  +  erf  [rj (7^2 ) 


<  0.0036( r 2  +  cr2). 


Therefore,  with  a  very  good  approximation,  r  «  <r>  and  cr2  «  <r2>  .  r 2 . 


To  compare  growth/coarsening  behavior  in  different  systems  and/or  at  different  conditions,  a 
nonnalized  particle  size,  p  =  r/r  ,  is  generally  used.  For  this  quantity,  the  PSD  function 
(Equation  16)  is  modified  to: 


(p-i):  1 

2  (cr/r)2y 


(17) 


Figure  1  shows  experimental  particle  size  distributions  for  a  number  of  alloy  systems.  These 
experimental  data  were  collected  in  ref.  [5]  and  are  shown  in  the  figure  as  open  circles.  The 
Llifshitz-Slyozov  distribution  function  is  shown  in  Figure  1  as  a  dotted  line.  It  is  highly  non- 
symmetric,  it  does  not  predict  particles  with  sizes  larger  than  1 .5  r  and  it  does  not  fit  well  the 
experimental  data.  At  the  same  time,  the  normal  distribution  function  describes  the  experimental 
data  rather  well  and  it  has  flexibility  to  change  the  PDS  width  and  height  by  assigning  different 
dr  ratios. 

Figure  2  shows  experimental  particle  size  distribution  for  Sn  particles  after  coarsening  for  40 
hours  in  a  Pb-Sn  liquid  [13].  Theoretical  PSDs  derived  from  different  modern  theories  of  the 
particle  coarsening  are  also  shown  in  this  figure.  These  PSDs  are  non-symmetrical  and  their 
agreement  with  the  experimental  PSD  is  rather  speculative.  The  blue  dashed  line  is  the  normal 
distribution  given  by  Equation  (17)  with  dr  =  0.4,  which  agrees  with  the  experimental  data 
very  well.  There  is  a  number  of  other  examples  showing  experimental  PSDs  developed  during 
particle  growth  and  coarsening  to  be  described  by  the  normal  distribution  function 
[14, 15, 16, 17, 18, 19,20,21],  in  agreement  with  our  analysis . 

In  conclusion,  our  theoretical  analysis  of  particle  growth  and  coarsening,  which  includes 
stochastic  flux,  predicts  nonnal  distribution  of  the  particle  sizes,  independent  on  the  particle 
volume  fraction  (see  Equations  16  and/or  17).  This  result  agrees  well  with  the  numerous 
experimental  data  and  it  solves  a  mysterious  problem  related  with  a  non-symmetrical  size 
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distribution  in  earlier  LSW  approach.  As  we  can  see  from  the  present  analysis,  taking  into 
account  a  stochastic  process  of  solute  diffusion  inside/from  the  particles  adds  an  additional 
(diffusion)  term  in  the  differential  Equation  (15)  for  the  PSD  function.  This  stochastic  process 
was  ignored  in  the  previous  LSW  theories,  which  led  to  an  artificial  requirement  that  the  LSW 
PSD  function  and  all  its  derivatives  should  go  to  zero  above  some  finite  value  of  the  particle 
size;  otherwise  the  integral  of  the  LSW  distribution  function  would  be  logarithmically  divergent 
[1,2,3]. 
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Figure  1.  Experimental  size  distributions  of  growing  particles  in  different  alloy  systems  (shown 
as  open  circles  and  taken  from  Figure  1  of  ref.  [5]),  Lifshitz  and  Slyozov  distribution  (dotted 
line)  [3]  and  normal  distributions  (solid  lines,  Equation  17)  at  a/<r>  =  0.2  (blue  line),  0.28  (black 

line),  and  0.4  (red  line). 
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Figure  2.  Experimental  (dark  solid  line)  and  theoretical  (other  gray  lines)  particle  size 
distributions  for  Sn  particles  by  Snyder  et  al.  [13].  The  theoretical  PSDs  indicated  in  the  legend 

are  from  modified  LSW  theories  and  the  blue  line 
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